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Abstract: In this paper we will present the methodology and pedagogy of Elementary 
Mathematical Modeling as a one-semester course in the liberal arts core. We will focus on the 
elementary models in finance and business. The main mathematical tools in this course are the 
difference equations and matrix algebra. We also integrate computer technology and cooperative 
learning into this inquiry-based learning course where students work in small groups on carefully 
designed activities and utilize available software to support problem solving and understanding of 
real life situations. We emphasize the use of graphical and numerical techniques, rather than 
theoretical techniques, to investigate and analyze the behavior of the solutions of the difference 
equations. 

As an illustration of our approach, we will show a nontraditional and efficient way of introducing 
models from finance and economics. We will also present an interesting model of supply and 
demand with a lag time, which is called the cobweb theorem in economics. We introduce a sample 
of a research project on a technique of removing chaotic behavior from a chaotic system. 

Key words: mathematical modeling, difference equations, cooperative learning, inquiry-based 
learning, cobweb theorem in economics. 


1. Introduction 

We introduce some models from finance and business represented by a single first order difference 
equation as a part of a course in mathematical modeling aimed at freshmen undergraduate students. 
The course is designed as a one-semester course in the liberal arts (general education) core, where 
students are required to take one or two mathematics courses. 

The main mathematical tools are difference equation, matrix algebra and Markov chains. The course 
does not require training in calculus. 

The use of technology in the form of computer algebra system (CAS) such as DERIVE and/or 
spreadsheet such as Excel, and cooperative learning [3] are integral part of the course. Students work 
in small groups on carefully designed activities that guide them in exploring and discovering the 
concepts on their own [4], It is important to note that most of the concepts are intuitively introduced. 

In this paper we show how the students use difference equations to model annuities, sinking funds, 
mortgages and amortization. They use a computer algebra system, graphing calculator, or spreadsheet 
to iterate the difference equations to find the solutions and investigate the behavior of the graphs. We 
will also present an interesting model of supply and demand with a lag time, which is called the 
cobweb theorem in economics. We introduce a sample of a research project on a technique of 
removing chaotic behavior from a chaotic system. 


This paper is an expanded version of the presentation given to the second Computer Algebra and Dynamic Geometry 
Systems in Mathematics Education (CADGME) conference at the University of Linz, Austria, in July 2009. 
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2. Some Models from Finance and Business 

2.1 Annuities. The students work in small groups of 4 on activities that are designed to allow them to 
model real life financial situations with difference equations and to explore related concepts on their 
own. 

Activity 2.1 It is known that, for certain income brackets, money deposited in an individual 
retirement Account (IRA) and its interest are tax-deferred. Karen Smith sets up an IRA at the 
beginning of the year and deposited $3,000. Karen will deposit $2,000 at the end of each year, and the 
interest rate is 6% compounded annually. 

A) Model this situation by a difference equation. 

B) What will be the future value of the account at Karen’s retirement age of 60 if she starts the IRA at 
the age of 25? 

Discussion: The students in groups work on modeling this situation and then use DERIVE or 
spreadsheet to iterate the difference equation to answer part B. After a few minutes the instructor asks 
the representative of one group to represent their approach and then allow discussion among groups. 
Here is a typical students approach. Let y n be the balance in the account after n years. This situation 
can be modeled by the difference equation: 

y n = T„_i + 0.06y„_, + 2,000 = 1.06> , n _ 1 + 2,000 ( 2 . 1 ) 

where y 0 = 3, 000. We are looking for y 35 , which can be obtained by iterating the difference 
equation (2.1). Most of the groups use DERIVE’s function ITERATE to find y 35 by authoring and 
approximating the following expression: 

ITERATE (1.06y + 2000, y, 3000, 35) 

The output is 245,927 which is y 35 . Note that the in order to obtain the ordered pairs («,»,), n = 0, 
1, 2, ..., 35, and graph y n vs. n, author, approximate, and plot the DERIVE expression: 

ITERATES ([n + 1, 1.06y + 2000], [n, y], [0, 3000], 35) 


Activity 2.2 Assume you open an account that pays 4% compounded annually and deposit $2,000. 
You will deposit 2% more into you account than you deposited in the previous year. 

A) Model this situation by a difference equation. 

B) What is the total amount in the account after 20 years? 

Discussion: The students in groups model this situation but the iteration of the difference equation 
is not a straight forward. The following is an approach by a group of students to model the situation. 
Let y n represent the total amount in the account after n years and let d n be the deposit at the end of 
the n th year. We have 


d n+1 =d n +0.02d n 


d n+l = 1 .02d n , </ 0 =2000 

(2.2) 

y„ + i =T n +°- 04 T„+< + i 


W+i = 1 -04W + d n+] , y 0 = 2000. 

(2.3) 


However, most of the students can not iterate the difference equation (2.3) since its right-hand side 
contains two dependent variables y n and d n+v The groups and instructor’s discussion suggests that 
equation (2.3) can be iterated if its right-hand side is transfonned into a function of the dependent 
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variable y n and the independent variable n. Students from previous sessions know that the difference 
equation x n = ax n _ x , where a is a constant, is equivalent to x n = a"x 0 . Consequently, the difference 
equation d n = 1 ,02d n _ l is replaced with d n = (1.02 )”d Q = 2000(1.02)". Thus, the current situation is 
modeled by the difference equation: 

y„ +l = 1 -04 y g + 2000(1.02)' ,+1 , = 2000 (2.4) 

Now DERIVE’s function ITERATE can be used to find the ordered pair (20,y 20 ) by authoring and 
approximating the expression: 

ITERATE ([n + 1, 1.04y + 2000(1.02) A (n + 1)], [n, y], [0, 2000], 20) 

The output is [20, 76310.10]. 

After these activities the instructor introduces the concept of annuity and the associated terminology 
such as ordinary annuity, the term and future value of an annuity. 

2.2 Mortgages and Amortization. While working on the following activity the students explore 
modeling loans and determine the periodic payments on a loan as well as other related questions. 

Activity 2.3 You want to buy a new car that costs $16,000. The dealer offers $4,000 for a trade-in 
with a down payment of $2,000. You borrow $10,000 from a bank for 5 years at an annual rate of 
interest 12% compounded monthly. The payments are made at the end of each month. What will the 
monthly payment be? 

Discussion: Students model the given situation by a difference equation such as: 

y„ =y„- 1 +0.01y n _ l -p = 1.01y n _ l -p, (2.5) 

where y n represents the balance of the loan after n months and p is the monthly payment, where the 
interest rate per month = 0.12/12 = 0.01. In this situation y 0 =10,000 and y 60 = 0. Since the value 
of p is unknown, we cannot iterate equation (2.5). The students usually estimate a value of p. If you 
take a loan of $10,000 without interest and you need to repay it over 60 months, the monthly payment 
would be 10,000/60 = $166.67. Since there is an interest rate on the loan, the monthly payment must 
be more than 166.67. Students simply select an estimated value of p, substitute it in (2.5), and use 

DERIVE ITERATE routine to evaluate y 60 , which should be 0 or very close to 0. The students adjust 
their estimation of phased on the value of y 60 . After several trials they find that if p = 222, 
y 60 = 36.29 and if p = 223, y 60 = —45.36. This concludes that $222 < p < $223. P ~ $222.25. 

Usually some students would wonder if there is a fonnula that can be used to solve the above problem 
instead of guessing the monthly payment. The students are guided to develop the analytical solution of 
a non-homogeneous first order linear difference equation in the fonn 

y„=ay n -i+b ( 2 . 6 ) 

where a and b are constant real numbers and a . We have y x = ay 0 +b, y 2 =ay l +b = 
a(ay 0 +b) + b = a 2 y 0 + ab + b . Similarly, y n = a"y 0 +a"~ l b + a"~ 2 b + ... + b = 
ci"y 0 + 6(1 + a + a 2 +... + a "~'). Using the sum of a geometric sequence, we get 

y„ = ~r + (To + “^rK if a * 1 ( 2 - 7 ) 

(7-1 (7-1 
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The instructor asks the students to use (2.7) to find the monthly payment p in the above activity. 
They substitute values of: a = 1.01, b = -p, y 0 = 10,000, n = 60 and y 60 = 0 in (2.7) and solve 
the equation for p, which is $222,444. 


2.3 Equilibrium and Stability. The following activity is designed to intuitively introduce students to 
the equilibrium value of a first-order linear difference equation with constant coefficients and the 
stability of the equilibrium value. 

Activity 2.4 For each of the following difference equations find the ordered pairs ( n,y n ), 
n = 0,1,2,..., 24 for the given initial conditions of y 0 and describe the graphs: 

A) y n = 0.6j„_j +8, ( n,y n ), n = 0,1,2,...,24 for (i) y 0 = 8, (ii) y 0 = 20, and (iii) y 0 = 28. 

B ) y„ = 1 - 1 T „-1 -4, (n,y„), n = 0,1,2,...,20 for(i) y 0 =35, (ii) y 0 = 40, and (iii) y 0 = 45. 

C) y„ = -0.6y n _ l + 16, (n,y n ), n = 0,1,2,...,24 for (i) y 0 = 6, (ii) y 0 =10, and (iii) y 0 = 14. 

D ) y n =-1.5j„_i + 5, ( n,y „), n = 0,1,2,...,12 for (i) y n =1, (ii) y 0 = 2, and (iii) y 0 = 3. 

Discussion: The students use DERIVE or a spreadsheet to explore the behavior of the numerical 
solutions of the given difference equations with given initial conditions. They will recognize the 
equilibrium value of the difference equation and determine whether the equilibrium value is stable or 
unstable. The instructor develops with students the definition of an equilibrium value of a difference 
equation and how to detennine it. The definitions of stable and unstable equilibrium values are 
developed too. In addition, the instructor discusses with students the equilibrium values in the context 
of financial and business models. 

Now the students are ready to do the following activity. 

Activity 2.5 Make a conjecture about the stability of the equilibrium value of the difference 
equation 

y n = a y„-i+ b 

in tenns of the constants a and b. 

Discussion: Lively discussion among the group mates and the groups will lead students to 
conclude that an equilibrium value is stable if — l<a<l, that is |a|<l. It is interesting that the 
stability of an equilibrium value does not depend on the value of b. 


3. A Model of Supply and Demand 

We will consider the fanning industry where a crop is produced once each year and the fanner decides 
on the amount of the crop to be planted based on the current year price of the crop. Note that the 
amount of crop planted this year is harvested next year and would represent the next year supply. Let 
S and D n be the number of units supplied and demanded in year n respectively and let P be the 

price per unit in year n . If the price of the crop is high in the nth year, the fanner anticipates that the 
price level will be maintained and will decide to plant a larger amount of the crop. The following year 
with a large harvest of the crop, the supply S n+1 exceeds the demand D n+X and will result in dropping 

the price of a crop P n+l in order to create more demand to equal the supply. As a result of the drop in 
the price P n+l the fanner will decide to plant a smaller amount of the crop. Consequently, the supply 
S „+2 be less than the demand D n+2 and will result in increase of the price P n+2 and the fanner 
will decide to plant a larger amount of the crop and this process repeats itself. 
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We will make the following assumptions: 

1. The supply of the crop S n+l in year n +1 is a linear function of the crop’s price in the 
previous year P n . The supply S n+l increases when the price P n increases. 

2. The demand of the crop D n in year n is a linear function of the crop’s price in the same year 
P n . The demand D n decreases when the price P n increases. 

3. In any year n , the market price P n of the crop is the price at which the demand D n equals the 
supply S n . 

The students use these assumptions to construct a model, investigate it, and make conjectures. From 
assumption 1 we have 

A + 1 = m s P n+C s (3.1) 

where m s and c s are positive constants. The constant m s is the slope of the line (11), which is one 
unit increase in price, produces an increase of m s units in supply. m s represents the sensitivity of 
suppliers to price. The constant c s is the S -intercept. From assumption 2, the demand equation is 

A ,=-m d P„+c d (3.2) 

where m d and c d are positive constants. Since m d is positive, the slope of the demand line is 
negative and means that a one unit increase in price produces m d units decrease in demand. m d 
represents the sensitivity of consumers to price. Assumption 3 asserts that 

D„=S„ (3.3) 

From (3.1), (3.2) and (3.3) we get 

~™d P n+ C d = m s P n-l +C S 

P n =aP n . 1 +b (3.4) 

where a = -m s / m d and b = (c d - c s ) / m d . 

The price at which the supply and demand lines intersect is called in economics the equilibrium price, 
that is the intersection of S n and D n . Equivalently the equilibrium price, call it P e , is the equilibrium 

value of difference equation (3.4). The equilibrium price is called stable if the prices converge to P e , 
otherwise is called unstable. 

Equation (3.4) is a first order linear difference equation with constant coefficients. The behavior of the 
solutions of (3.4) depends on the value of the constant a. Noting that (7 < 0 , we are interested in 
investigating (3.4) in three cases: Case 1: -1 < a < 0 (i.e. m s < m d ); Case 2: a — -1 (i.e. m s = m d ); 
Case 3: a <-1 (i.e. m s > m d ). 

Activity 3.1: Consider the supply and demand model (3.4) in the following cases: 

Case 1: m s = 0.25, m d = 0.5, c s =2, c d = 8, and .P 0 = 10 

Case 2: m s = 0.5, m d = 0.5, c s =1, c d = 8, and P 0 = 8 
Case 3: m s = 0.6, m d = 0.5, c s =1, c d = 7.6, and P 0 = 6.5 
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For each case do the following: (i) Find the equilibrium price P e . Find a numerical solution 
(n,P i: ),n = 0,1)•••,14, graph it and analyze its behavior, (ii) Detennine whether the equilibrium price 
P e is stable. Describe the relationship between stability/instability and sensitivity of the suppliers and 
consumers to price changes. 

Discussion: In Case 1, a = -m s l m d =- 0.5, b = (c d -c s )/m d = 12, and equation (3.4) 
becomes, P n+1 =-0, 5P n +l2. The equilibrium price P e =b/{\ -a) = 8. The graph of the ordered 
pairs ( n,P n ),n = 0,1,...,14 is shown in Figure 3.1. The price oscillates around 8, that is the solution is 
alternately above and below the line P = 8. The amplitude of the oscillation decreases and eventually 
the price P n approaches 8 as n becomes large. Consequently the equilibrium price is stable. 

In Case 2, a = -1 ,b = 12, and P e = 6. The graph of ( n,P n ),n = 0,1,...,14 is shown in Figure 3.2, 
where the price oscillates between 8 and 4 that is between P 0 and —P 0 4- b . 

In Case 3, a = -\.2,b - 13.2, andP e =6. The graph of (n,P n ),n = 0,1,...,14 is shown in Figure 3.3, 

where the price oscillates around P = 6. The amplitude of the oscillation increases without bound. 
The price diverges from the equilibrium price. Consequently, the market price is unstable. Note that 
this model will fail if the price becomes negative. 
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Conjecture: Suppose that the price of a product satisfies the above three assumptions in the supply 
and demand model. It can be concluded from the above activity that the stability of the equilibrium 
price, that is the market price, depends on the value of a = —m s / m d . We have three situations: 

Situation 1: If o >-l, that is m s < m d , then we have a stable price, i.e. the market price tends to 

stabilize. In other words, we have a stable market if the sensitivity of suppliers to the price is less than 
the sensitivity of consumers to the price. 

Situation 2: If a < — 1, that is m s > m d , then we have an unstable equilibrium price. This means 

that we have an unstable market if the sensitivity of suppliers to price is greater than the sensitivity of 
consumers to the price. 

Situation 3: If a — — 1, that is m s = m d , then the market price oscillates and is a 2-cycle. The 
market price is unstable. 

These three cases can be summarized in the following statement: If the producers are less sensitive to 
the price than the consumers, then the market price tends to stabilize. Otherwise, we have an unstable 
market price. 

Practical considerations: In practice, it is difficult to calculate the constants in this model, 
m s ,m d ,c s , and c d for a particular product. In reality, we are interested in knowing whether the 
equilibrium price is stable or not. Recall that the above conjecture concludes that the equilibrium price 
is stable if m s <m d . Consequently, to determine the stability of the equilibrium price, we need to 

compare m s and m d . In other words, we need to observe the reaction of suppliers and consumers to 
different prices of the product to determine who the more sensitive one is. 

A big concern arises if it has been determined that n? s > m d for a certain product, for example 
potatoes, and consequently, the market price is unstable. In order for the market to change and become 
stable, m s < m d is needed. One way to increase m d (the consumers’ sensitivity to price) is for the 

government to act as the consumer. This can be done in one of two ways: (i) the government supports 
the price; (ii) the government pays the fanners to limit the potatoes they grow. 


4. Logistic equations and chaos. 

Many situations cannot be modeled by linear first order difference equations in the fonn (2.2), but can 
be modeled by first order nonlinear difference equations in the fonn (4.1), 

x„ = x n _ ] +ax n _ ] -bx 1 n _ ] (4.1) 


where a and b are constants. An equation in the form (4.1) is called logistic difference equation. 
1 -(- ci 

Letting x n =- y n , equation (4.1) is scaled to equation (4.2) which is called logistic map. 


T„=(l + «)T„-i(l-T„-,) (4-2) 

Students explore the behavior of the solutions of (4.2) for different values of a,0 < a <3. They 
recognize the intervals for a when solutions behave in different ways such as: the solution oscillates 
around the limiting value and converges to it (1 < a < 2). This type of oscillation is called damped 
oscillation ; the solution oscillates between two values (2 < a < 2.4495...) and is called periodic 
solution with period 2 (2-cycle)', the solution repeats itself every four iterations 

(2.4495... < a < 2.56...) and is called periodic solution with period 4 (4-cycle)', and the solution does 
not approach a fixed point or a periodic cycle and behaves erratically and is very sensitive to initial 
conditions (2.56994... < a < 3). This type of behavior is called chaotic. 
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5. Sample research project: a technique to remove chaotic behavior from a system. 

There are some methods to control chaotic systems. Giiemez and Matias [2] suggested an interesting 
and simple method to control detenninistic chaos of a system represented by a single difference 
equation. The method works by applying a series of proportional feedbacks (adjustments) on the 
system variable performed periodically with a certain period. The method, like other methods, 
stabilizes a given unstable orbit. However the method is simpler than other methods. 

Consider the logistic map (4.2), with A = 1 + a , 

y n (5.1) 

The control algorithm consists of the application every A n iterations of a feedback to the variable y 
having the fonn 

y n = y,M+y) (5-2) 

where y is a real number that represents the strength of the feedback, and — 1 < y < 1. 
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Figure 5.4 The control algorithm is applied to (5.1) with A = 3.9,/ = —0.2, and An = 4. 


Let An be some positive integer representing the period of control. Thus the control is given by, 

y„ + 1 = Ay n (1 - y n ) if n does not divide An, otherwise y n+l = Ay' n (l~y„) (5.3) 

One way to iterate the above function (5.3 ),y n+l , is to define it as piecewise function. We will use 
DERIVE function CHI. Recall that CHI (a, x, b) = 1, if a < x < b. If x < a < b or a < b < x, CHI (a, x, 
b) = 0. We provide two examples of applying this control algorithm to different values of A, An, and 
V- 

Example 1. Figure 5.1 shows the graph of y n versus n for the logistic map (5.1) with A = 3.7, 
while Figure 5.2 shows the graph of (5.1) after applying the control algorithm with y = 0.05, and 
An = 5. The ordered pairs ( n,y n ),n = 1,2,...200 may be obtained by the following expression in 
DERIVE: 

ITERATES ([n+1, 3.7y(l -y) CHI (0.1, MOD (n, 5), 4.1) + 3.7y(l + 0.05) 

(1 -y(l + 0.05)) CHI (-0.1, MOD (n, 5), 0.1)], [n, y], [1, 0.4], 200) 

Example 2. Figure 5.3 shows the graph of versus n for the logistic map (5.1) with A = 3.9. 
Figure 5.4 shows the graph of (5.1) after applying the control algorithm with A — 3.9, y = —0.2, and 
An = 4. The ordered pairs (n, y n ), n = 1,2,...650 is obtained by the following DERIVE expression: 

ITERATES ([n+1, 3.9y(l -y) CHI (0.1, MOD (n, 4), 3.1) + 3.9y(l - 0.2) 

(1 -y(l - 0.2)) CHI (-0.1, MOD (n, 4), 0.1)], [n, y], [1, 0.2], 650) 
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Figures 5.2 and 5.4 show that the periodicity of a system does not set immediately after applying the 
control algorithm. Note, that this control algorithm does not change the system parameters which is an 
advantage. However, the control algorithm does not work for every values of y and An, which is a 
disadvantage. 

Giiemez and Matias algorithm can be applied to several real unstable systems such as the economy 
which is a complex system [1], In case of economy, the proportional feedbacks in this control 
algorithm can be represented by injection, if y is positive, or withdrawal, if y is negative, of 
macroeconomical quantity. 

The student research teams are required to review Giiemez & Matias’ and other related articles. They 
find values of y and An where the control algorithm is valid and determine the algorithm’s 
limitations. Each team researches possible application of this algorithm in a specific area and presents 
their findings to the rest of the teams. 

Conclusion 

Modeling real life situations excites students and shows them the relevancy of abstract mathematical 
concepts to their lives, and helps convey the elegance and wide applicability of mathematics. 
Difference equations represent a very sophisticated and powerful mathematical tool to model a wide 
range of real life discrete time situations in diverse areas. Moreover, this powerful tool does not 
require sophisticated mathematics background, being accessible to anyone who has successfully 
completed high school algebra. 

Research in mathematics education shows that students learn mathematical concepts more 
constructively by making mental images. Easy to learn and use symbolic, numeric, and graphic 
software helps students in constructing these mental images in a creative and productive way. 
Moreover, the use of computers frees students from tedious calculations and associated with these 
calculations boredom, which in turn allows students to focus on translating a problem into 
mathematical notations, finding a solution, interpreting the numerical and the graphical infonnation 
provided, and then making conjectures and writing about their findings and observations. With the use 
of computers, real life and realistic problems can be considered in class/lab activities. 
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